Environment Assisted Metrology with Spin Qubits 
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We investigate the sensitivity of a recently proposed method for precision measurement [Phys. 
Rev. Lett. 106, 140502 (2011)], focusing on an implementation based on solid-state spin systems. 
The scheme amplifies a quantum sensor response to weak external fields by exploiting its coupling 
to spin impurities in the environment. We analyze the limits to the sensitivity due to decoherence 
and propose dynamical decoupling schemes to increase the spin coherence time. The sensitivity is 
also limited by the environment spin polarization; therefore we discuss strategies to polarize the 
environment spins and present a method to extend the scheme to the case of zero polarization. The 
coherence time and polarization determine a figure of merit for the environment's ability to enhance 
the sensitivity compared to echo-based sensing schemes. This figure of merit can be used to engineer 
optimized samples for high-sensitivity nanoscale magnetic sensing, such as diamond nanocrystals 
with controlled impurity density. 



I. INTRODUCTION 

Quantum metrology seeks to achieve precision mea- 
surements with an accuracy beyond the limits imposed 
by the central limit theorem [1] (the standard quantum 
limit, SQL). Although many proposals for achieving the 
quantum limits of sensitivity (as defined by the Hciscn- 
berg bounds) have been presented, they arc often diffi- 
cult to implement in practice. The main challenges arise 
from the deleterious effects of noise and decoherence on 
the (entangled) states required for quantum metrology 
and from the unavailability of the Hamiltonians and mea- 
surement strategies needed to create and readout these 
entangled states. 

We recently introduced a scheme [2] that aims at over- 
coming these two challenges. We proposed to use the 
environment of the sensor as an additional resource for 
metrology and we showed how to achieve the desired in- 
teraction Hamiltonian using coherent control techniques. 
In this paper we focus on one possible implementation of 
this environment assisted metrology (EAM) scheme - a 
spin sensor embedded in a bath of other spins - in order 
to derive more detailed results on the sensitivity achiev- 
able. In addition, we will analyze in depth the effects of 
decoherence and of finite polarization. 

The paper is organized as follows. In Section II we 
present the EAM scheme: the control sequence that 
achieves it and the sensitivity gain in the idealized sit- 
uation of no decoherence. This restriction is lifted in 
Section III, where we analyze the effects of decoherence, 
both analytically and with numerical simulations. We 
further provide strategies to reduce the effects of decoher- 
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ence. In section IV we use these results to derive limits of 
the proposed EAM strategy and compare them to usual 
strategies that do not take advantage of the environment. 
Since the sensitivity depends on the polarization of the 
spin environment, we propose in section V schemes for 
polarizing these ancillary spins and we further extend 
the scheme to the case where no polarization is available. 
A second extension of the EAM method is presented in 
section VI, where more general spin systems are studied. 



II. THE ENVIRONMENT ASSISTED 
METROLOGY SCHEME 

We consider the metrology task of measuring a param- 
eter b via its interaction with a quantum probe. The task 
can be accomplished by using a Ramsey scheme, where 
a two level system is first prepared in a superposition 
of the two states, which then acquire a phase difference 
that is mapped onto the populations by a second pulse. 
An example of this scheme is magnctomctry with solid- 
state spins [3], where the probe interacts with the exter- 
nal magnetic field via a Hamiltonian rl oc bS z , acquiring 
a phase cx bt during the interrogation time t. Then, the 
bound to the sensitivity is set by the dephasing rate that 
limits the time the probe can interact with the external 
field associated to the parameter to be measured. 

Coherent control techniques can be used to isolate the 
probe from its environment, thus increasing the coher- 
ence time. If the environment interacts as well with the 
external field to be measured - as it is the case for a spin 
bath - a different strategy is possible: in Rcf. [2] we 
showed that the spin environment can be used as a re- 
source in this case, by mapping the phase acquired by the 
environment spins onto the probe spin before readout. 
Here we provide more details of the method presented in 
Ref. [2] and consider several extension of the work. To 
this end we assume that the spin environment can be col- 
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lectively controlled and partially polarized. These spins 
could thus be considered as an ancillary system. Still, 
since they cannot be addressed individually nor read out, 
they cannot be used directly as probes or in sequential 
adaptive schemes [4-6]. In addition, because their cou- 
plings to the probe spin cannot be switched off, they are 
a cause of decoherence for the probe spin (as we will see 
in section III) and thus they can be considered as envi- 
ronment. Nevertheless we show that one can make active 
use of these spins, to increase the sensitivity of a mea- 
surement. 

Ancillary qubits have been considered as a resource for 
parameter estimation [7] in a scheme inspired by the de- 
terministic quantum computing with only one pure qubit 
(DQC1) model [8]. In that scheme, the probe qubit is ini- 
tially prepared in a superposition state, then the ancillary 
system interacts with the external parameter conditional 
on the state of the probe, which is finally readout (see 
Fig. l.a). When the conditional evolution is given by 
the operator U = e~ lbt l * (where P* are the ancilla 
spin operators) the sensitivity achieves the SQL (scaling 
as 1/ y/n where n is the number of ancillary qubits) for 
ancillas in a completely mixed state [7] and the Heisen- 
berg limit for pure state (scaling as 1/n). In that case, it 
is convenient to read out the y-componcnt of the probe 
spin, which gives a signal S = sin(n&£). Since the signal 
is enhanced by a factor of n for small fields nbt <C 1 this 
yields an Hciscnbcrg-limitcd sensitivity scaling as 1/n. 
Indeed, for pure input states, the circuit creates an en- 
tangled state that provides a signal enhancement. Below 
we modify this scheme so that it can be implemented for 
realistic physical systems. 

In general, the ancillas dependence on the external pa- 
rameter cannot be controlled by the probe spin, as it is 
implicitly assumed above. Thus it is necessary to inter- 
sperse the evolution under the interaction with the exter- 
nal field with C-Not gates (Fig. l.b). With this modifi- 
cation we achieve a similar evolution as before. However, 
even this simpler scheme cannot be easily implemented 
and is not compatible with our assumptions of limited 
control on the environment spins: if the ancillas arc spins 
in the environment, it is not possible to control them in- 
dividually, thus the C-Not gates cannot be implemented 
since the required interaction time for the C-Not oper- 
ations will be different for the different spins. The key 
to using the environment spins - with the corresponding 
limited control - as a resource for parameter estimation 
is to realize that the scheme works also if the the con- 
trolled gates are not ideal 7r-rotations. The rotations can 
differ for different spins, as long as the state of the probe 
spin is flipped (Not gate) before the second set of con- 
trolled gates: this ensures that all the environment spins 
contributes constructively to the final phase, as we derive 
below. 

We note that the spin flip of the probe achieves two 
other results: first, it makes the evolution insensitive to 
static noise (as produced for example by a very slowly 
varying spin bath) since the gate amount to a spin echo 
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FIG. 1: (a) Ideal circuit for EAM with ancillary qubits, based 
on a DQC1 scheme, (b) More realistic circuit, where the in- 
teraction with the field to be measured is not conditional on 
the probe spin. The controlled- I x gates denote C-Not gates, 
in the ideal model, and are reduced to more general rota- 
tions (with different angles for different spins) in the realistic 
scheme. 



for the probe spin. Secondly, the echo pulse refocuscs 
the entanglement created in the first half of the circuit; 
this operation cancels undesired terms in the signal that 
would arise when considering a more realistic scenario 
where both the external field and the couplings to the 
probe spin used to create controlled rotations are always 
present at the same time. 

The idealized scheme in Fig. 1 can be implemented 
in practice with realistic resources, with the EAM pulse 
sequence of Fig. 2. We consider a system comprising a 
sensor spin (S = 1) and environment spins (/ fe ), which in 
a convenient rotating frame on resonance with the m s = 
0—>l transition, is described by the Hamiltonian: 



n = b(t) ( ls s z + 7/ E fe + E A fc s,J* 



|0)(0| 



L&(*)7/E*^] + 
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(1) 



where b(t) is the external field to be measured, "fs.i arc 
the gyromagnetic ratios of the probe and environment 
spins respectively, Xk are the dipole couplings between 
the sensor and environment spins, and |0) (|1)) denotes 
the m s = (m s = 1) eigenstate of the S z operator. 

We choose a spin-1 system for its analogy with 
Nitrogen- Vacancy centers in diamond [9, 10] as they 
have emerged as good quantum probes of magnetic 
fields [3, 11, 12] for their controllability, optical readout 
and long coherence times. In addition Nitrogen paramag- 
netic impurities (PI centers [13]) can act as the environ- 
ment spins, since they can be collectively controlled [14]. 
The choice of a spin-1 system is in addition important 
since the presence of an eigenstate with zero eigenvalue 
effectively allows shutting off the interaction between the 
probe spin and the environment spins at given times: this 
flexibility makes the EAM scheme easier to implement. 
We will lift this restriction and examine more general 
case in section VI. 

In the sequence in Fig. 2 the probe spin undergoes a 
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FIG. 2: EAM pulse sequence: the vertical bars represent mi- 
crowave pulses on resonance with the probe (top part of the 
figure) or environment spins (center), performing the labeled 
rotations. We assume that the field to be measured is an AC 
field synchronized with the pulse sequence as shown in the 
bottom of the figure. 



spin-echo sequence induced by pulses on resonance with 
the transition between the states |0) and |1) before being 
measured. For any given evolution of the environment, 
the signal can then be calculated from S{t) = [1 +£>(£)] /2, 
with [15, 16]: 

S{t) = Im [lr [u U lPcm Ulul}~\ (2) 

Here the propagators Ui = e~ tUit are defined as the evo- 
lution of the environment spins in the m s = i manifold, 
where H = b(t)7i J2 k 4* and Ui = 7s&(*)+Ek(7z&(*) + 
Xk)Iz ( see Eq. 1). The pulsed evolution of the environ- 
ment, giving the propagators Ui, can be most easily cal- 
culated in the toggling frame [17], the interaction frame 
defined by the control pulses. In this frame, the Hamil- 
tonian (1) becomes piecewise time-dependent, with op- 
erators alternating between the z- and x-axis. 

The evolution for the sequence of Fig. 2 and the result- 
ing signal (Eq. 2) can be calculated exactly in the case 
of a single ancilla. Here we will present only the result 
for many ancillas in the limit of small field 6, following 
the derivation of Rcf. [2]. We neglect for the moment 
the coupling of the sensor spin to the magnetic field and 
only keep first order term in the field. By expanding the 
exponentials, the only terms contributing to the signal 
arc then 
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inserted dt in integrals where B2 = 
The signal is then given by S = \ [1 



$ = ls B l 



1 



- sin($)], with 
8 



(3) 



where B\ = i I J Q 2 b(t)dt — J r b(t)dtj is the contribution 

from the direct coupling of the sensor with the field and 
we have introduced the polarization P < 1 of the envi- 
ronment spins, so that the initial state of each spin in 
the environment is p k = 1/2 + PI^ . The factor in the 
square bracket is the amplification attained as compared 
to magnetometry performed via a spin echo [3]. We can 

always get an amplification, as sin (§-) is non- negative 
and changing the pulse phases always ensures that 7/-B2 
and 7s B\ have the same sign. 

For values of the couplings such that |Afcr| > 27r, or 

strongly coupled environment spins, the terms sin (-g-) 
average to \ . Weakly coupled environment spins (A& < 1) 
contribute instead with a factor tx and we obtain a 
total phase 

* = 1s~B~It 1 + P^=l 

isBi 



2^'(A fe r/8) : 



, (4) 



where n sc is the number of strongly coupled spins and 
the primed sum runs only over the weakly coupled spins 
(this last term can generally be neglected compared to 
the strongly coupled spin contribution). 

The sensitivity of the EAM scheme is easily calculated 
by noting that ideally the only noise contribution is the 
shot noise of the spin probe. For 75 = 7/ = 7 and assum- 
ing an oscillating field in phase with the echo sequence 
b(t) = bo sm(2irt / t) , the sensitivity [19, 20] per unit time 



C7(2+iPn sc )^ 1 ' 



(5) 



where we introduced the factor C [3] to include any non- 
ideality of the measurement procedure (here we assumed 
T = Nt, with N the number of repetitions of the mea- 
surement). The sensitivity scales as l/n sc achieving a 
Heisenberg-like scaling 1 . 

We note that even in this ideal case, there are two fac- 
tors that reduce the sensitivity: a limited polarization 
of the environment spins and the reduction of the time 
during which the interaction with the external field is 
effective (because of the scheme proposed, a phase is ac- 
quired which is proportional to only l/4 t/l of the total 
sequence time). 

The EAM scheme thus demonstrates that it is pos- 
sible to attain nearly Heisenbcrg limited sensitivity for 
metrology with a new class of entangled states (other 
than squeezed or GHZ states) that as we will see in the 
following are more robust to decoherence. Furthermore, 
these states can be created with limited control resources, 
thus opening the possibility of using spins in the environ- 
ment as a resource for metrology. 



1 Equation (5) is valid only for P ^ 0. We will analyze the case 
P = in Section V. 
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FIG. 3: Left: Simulations of signal decay for spin echo sequence (Red, dotted and dash-dotted lines) and EAM sequence 
(Black, dashed and solid lines). Right: Simulations of the two sequences with a WAHUHA sequence embedded in each time 
period (for 1 to 50 cycles). The time was normalized by the largest sensor-environment spin coupling, r ~ [7r/A ma x]- 
The dotted and dashed lines correspond to a 6% spin density and the dash-dotted and solid lines to a density « 1/8. In the 
first case about 25 environment spins were placed around the probe spin on a diamond lattice, while in the second case, about 
50 spins were simulated. The polarization of each environment spin was P = | . We took an average of 100 spin distributions 
to obtain mean decay values and performed each simulation using the disjoint cluster method [18], with clusters of 6 spins. 



III. DECOHERENCE 

The results in the previous section did not take into 
account the effects of decoherence caused both by the en- 
vironment spins used as an ancillary system and by any 
other residual bath. In this section we will take these 
effects into account and show that even in the non-ideal 
case the EAM sequence can provide a sensitivity enhance- 
ment with respect to other control scenarios (such as a 
spin-echo) that only aim at refocusing the interaction of 
the probe spin with the environment spins. 

A. Decoherence induced by the environment spins 
couplings 

The interactions among environment spins hamper the 
EAM scheme in two ways. First, flip- flops of environment 
spins lead to a loss of coherence of the probe spin. This 
effect is the same that is observed during a spin echo, 
and we will show that the resulting coherence time Ti 
is on the same order for the two sequences. Second, the 
interactions will also cause the environment spins to lose 
their internal phase coherence, resulting in a smaller ac- 
cumulated phase $. Still, this effect happens on a time 
scale t i given by the environment spin correlation time, 
which is usually longer than the probe coherence time, 
77 > Ti. Thus the sensitivity is ultimately limited by T2, 
as in the spin-echo case. 

Consider the system evolution as given by Eq. 2 (for 
simplicity in the absence of the magnetic field b). Now 
the propagators are given by the Hamiltonian 

H = b(t) ( ls s z + 7/ E 4 fe ) + E *kSj k ; 

where f^ij fixe the intra-bath couplings given by the mag- 
netic dipole interaction among spins. Because of the pres- 



ence of the couplings, the evolution in the two halves of 
the sequence is no longer the same, thus the interaction 
between the probe spin and the environment spins can no 
longer be perfectly refocused. This effect, usually called 
spectral diffusion, is observed as well in spin echo experi- 
ments and lead to the coherence time T2 . The addition of 
a modulation of the environment spins is not expected to 
change substantially the coherence time, as hinted by the 
short time evolution expansion presented in Ref . [2] . An 
exception is for a perfectly polarized bath: in that case, 
flip-flops are quenched in the spin-echo, but they are still 
allowed in the EAM scheme since they are enabled by 
the rotation of the spins during the protocol; the effect 
of flip-flop quenching is however noticeable only for very 
high polarization of the bath [21, 22]. 

From this argument we expect that one can have a sim- 
ilar interrogation time r in Eq. (5) for the EAM scheme 
considered here as for a simple spin echo sequence. Un- 
like for different entangled states [23], the enhancement 
from entanglement is therefore not counterbalanced by 
a decrease in the interrogation time r, and the EAM 
scheme does allow for a significant improvement of the 
sensitivity. 

We further verify this claim by simulations. We used 
the disjoint cluster approximation [18] to simulate the se- 
quence in Fig. 2 for a system comprising the probe spin 
surrounded by an environment of 25-50 spins randomly 
positioned in a cube with sides of unit length. By aver- 
aging over many spatial distributions of the environment 
spins, the simulation converges quickly even for small 
cluster sizes and it gives information about the average 
coherence time [24, 25]. 

The system we consider is inspired by a NV center in a 
nano-crystal of diamond in the presence of PI Nitrogen 
impurities [2], but the results are more generally valid. 
For comparison, we also simulated the evolution under 
a spin echo sequence. From the results in Fig. 3 we see 
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that the coherence time is not qualitatively different for 
the two sequences. The figure shows in addition that the 
coherence time depends on the density of the environ- 
ment spins, a fact that will be important in evaluating 
the sensitivity achievable with the EAM scheme. 

The second effect of the intra-bath couplings is to make 
the environment spin themselves loose their coherence, in 
a time on the order of their correlation time 17 , which is 
given by the rate of spin flip-flop driven by the dipolar 
interaction. If the environment spins are no longer in 
a coherent state, the phase they acquire does not add 
up constructively, resulting in a smaller phase <£. Still, 
this effect is comparable to the previous one, since the 
correlation time is at least on the same order of T2. 

In addition to the environment spins that are used as 
ancillary sensors, the system could be in contact with an 
additional spin bath. For example, in the case of the NV 
center in diamond this bath is given by the 13 C nuclear 
spins. The effects of this quasi-static bath are refocused 
by the tt pulse on the NV center and by the two tt/2 
pulses on the environment spins, which amount to a so- 
called "Hahn echo" [26] sequence. Any residual decay is 
again comparable to what is observed in a simple spin 
echo for the probe spin. 



B. Dynamical Decoupling 

An increase in the effective correlation time of the en- 
vironment spins would be beneficial in two ways, by both 
increasing the coherence time of the probe spin, through 
its influence on the sensor spin T2-time, and directly by 
improving the environment spin coherence. Dynamical 
decoupling schemes could achieve this goal. The dipolar 



Probe 
spin 

Environment 
spins 

Magnetic 
field 

WaHuHa 



T/2 



HI 

2 be 



1 

K \x 

1 



T/2 



hi 

2\y 





1 2iy 




it i 

2 \y 







T/2---- 






II 


2t 


II 


H= 


— JT 1 Jtl 

2ix 2\y 

z y 


X 


tci —jti 

z\y 2\x 
-y -z 



approximated by its average. A simple symmetrization 
of the pulse sequence [28] can further cancel out the first 
order correction, leaving errors that are only quadratic 
in the product nt c (and do not depend on the total evo- 
lution time, that could be given by many cycles) [17]. 

Fig. 4 shows how to incorporate a WAHUHA sequence 
within the EAM sequence. We modified the phases of the 
pulses with respect to the original sequence in order to 
obtain an effective coupling between the probe and envi- 
ronment spins oc -^S z ^2 Afe/^-j in the odd(even) time 
intervals. These phase changes do not affect the average 
of the dipolar Hamiltonian and hence the performance 
of the WAHUHA sequence. Unfortunately the modula- 
tion does not only averages out the dipolar Hamiltonian, 
but it also reduces the linear terms, by a factor l/v3- 
In many cases, the increase in coherence time more than 
compensate for this weighting factor. In Fig. 3 we sim- 
ulated via the disjoint cluster method the coherence of 
the EAM and spin-echo sequences, while applying the 
WAHUHA sequence in between the pulses. Comparing 
the results obtained in the absence of dynamical decou- 
pling, we see that the sequence is very effective in in- 
creasing the coherence time. 

A different strategy for directly increasing the probe 
spin T2 is to use more than one 7r-pulse during the total 
sequence time [3] . This technique is inspired by concate- 
nated dynamical decoupling schemes and in particular by 
the CPMG sequence [29, 30]. More generally, these ex- 
amples indicate that the EAM scheme can be combined 
with various forms of decoupling. 



IV. SENSITIVITY 

In the previous section we saw that the coherence time 
(and hence the time during which the phase can be ac- 
quired) depends on the density of the environment spins. 
For the EAM sequence, the signal too depends on the 
environment spin density, since it determines how many 
environment spins are close enough to the probe spin 
to be considered "strongly coupled". Thus the optimal 
sensitivity arises from a compromise between the envi- 
ronment spin density and the interrogation time. Includ- 
ing the probe decoherence due to the environment spins, 
as well as other bath contributions, yielding a coherence 
time Trf , the sensitivity of Eq. 5 becomes: 



FIG. 4: (color online) Embedding of a WAHUHA sequence 
in the EAM sequence. The WAHUHA is shown at the bot- 
tom, together with the "direction" of the Hamiltonian in the 
toggling frame. 

Hamiltonian can be refocused using homonuclear decou- 
pling sequences such as the WAHUHA sequence [27] . The 
pulse modulation gives a time-dependent Hamiltonian for 
the spin-spin interaction that averages to zero over a cy- 
cle time t c . If the modulation is fast compared to the 
couplings, the effective Hamiltonian over the cycle is well 



7re (-r/T 2 ) 3 e (T/T 2 B ) 3 

C 7V ^(2 + \Pn sc ) 



(7) 



The functional form we assumed for the decay is inspired 
by the measured behavior of NV centers in diamond [10, 
14] and usually arises from a Lorentzian spectrum of the 
bath. 

The couplings between the probe and environment 
spins scales as Aj, ~ l 2 l T \ (assuming dipolar interac- 
tion and setting 7/ = 75 = 7 for simplicity), with the 
distance to the probe spin. Then, for a fixed duration r 
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of the EAM sequence, the number of "strongly coupled" 
spins n sc scales as n sc {r) ~ 7 2 /ot, where p is the density 
of the environment spins. The probe coherence time also 
scales with the density as X2 oc 1/p. 

The sensitivity is then a function of two parameters: 
how many polarized spins are strongly coupled in the 
coherence time T2 and how much the coherence time is 
reduced with respect to the background bath coherence 
time by introducing the ancillary environment spins. We 
define a quantity Q — Pp-f 2 T2, which describes the "qual- 
ity" of the environment spins. A second quantity describ- 
ing the reduction in coherence time due to the ancillary 
environment spins is given by the ratio r = T 2 B /T2- 

The EAM sensitivity then depends only on these two 
parameters and the bath coherence time, such that 
Eq. (7) becomes 



ne (l+r 3 )(r/T 2 B ) 3 

C 7 ^(2 + £*rQ) 



(8) 



We can further optimize the sensitivity with respect to 
the interrogation time r and compare it to the case where 
the field is measured by a probe spin (via a spin-echo 
sequence) in the presence of the background spin bath 
only (that is, no environment ancillary spins). In this 
case, the sensitivity is given by [3] 



^?echo,l 



ire 



(r/T 2 B f 



(9) 



As show in Fig. 5, the EAM sensitivity as given by 
Eq. (8) improves up to r — 1, where the decoherence 
due to the environment spins becomes more important 
than the background bath. The improvement depends 
on the "quality" Q, since for higher Q there are more 



strongly coupled spins at a given T2 time. It is then 
clear that there is an optimum number of environment 
spins that one would want to introduce in the system to 
obtain the optimal sensitivity. Alternatively, the quality 
Q can be improved by increasing T2 using the dynamical 
decoupling methods we introduced in section IIIB. 

If the number of environment spins is instead fixed, 
we are interested in comparing the EAM and spin-echo 
scheme for a given system (e.g. a given nanocrystal of 
diamond). In Fig. 6 we plot the ratio of the EAM sensi- 
tivity to the spin-echo sensitivity : 



TIC 



(l+r 3 )(r/T 2 B ) 3 



2C 7v ^ 



(10) 



(in the figure this expression is optimized with respect to 
the interrogation time r). In the high r limit the sensi- 
tivity ratio depends only on the Q factor, as 77/^echo ~ 

(1+2-V3Q) • 

To estimate the potential sensitivity improvement of 
the EAM method we express Q in terms of measurable 
quantities. Specifically, we can write the environment 
quality as Q = Ppj 2 T 2 w P\/J2 x l T 2, where we used 
the fact that 7 2 /3T = ti sc (t) 
distribution of couplings M2 



The average 
A 2 , is related to the 
second moment of the probe spin, which gives its dephas- 
ing time AI2 = l/TJ. Then the sensitivity improvement 
is given by the ratio T2/T2, which can be quite large in 
many system. 



V. POLARIZATION AND SENSITIVITY 

The sensitivity discussed in the previous section de- 
pends on the polarization of the environment spins. In 




FIG. 5: (Color online) Sensitivity of the EAM scheme nor- 
malized by the sensitivity of a spin-echo scheme in the absence 
of any environment spin, rj/^echo.l- The curves corresponds to 
quality factors Q = 10 (dotted), 20 (solid), 30 (dashed) and 
50 (dash-dotted). The ratio improves until r = /T2 = 1, 
where the decoherence induced by the added environment 
spins overtakes the background decoherence. 
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FIG. 6: (Color online) Sensitivity ratio rj/rj BC ho between EAM 
and spin echo schemes as a function of r = T 2 S /T2 for quality 
factor Q = 10 (dotted), 20 (solid), 30 (dashed), 50 (dash- 
dotted). When the decoherence induced by the added envi- 
ronment spins dominates (large r) the sensitivity ratio is only 
determined by the environment quality Q. 
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this section we first propose methods for creating this po- 
larization, under the assumption that the probe spin can 
be polarized at will. We then generalize the EAM scheme 
to the case where no polarization is available. This gener- 
alization will furthermore prove useful in the case where 
the field to be measured is affected by a random phase. 



A. Polarizing the environment spins 

In an environment assisted magnetometer working at 
room temperature, the environment spins will be in a 
thermal state, close to the maximally mixed state. Polar- 
ization need then to be created by relying on the probe 
spin and the Hamiltonian (Eq. 6) that is required for 
the measurement scheme. To do this we assume that 
the probe spin can be repetitively polarized: this is the 
case e.g. for an NV center that can be polarized op- 
tically. Polarization could then be transferred to the 
spins in the environment by a swapping Hamiltonian such 
as Hsw ~ (S x Ix + Syly). Although this operator is 
contained in the dipole-dipole Hamiltonian, it is usually 
quenched in the rotating frame, if the energies of the two 
spin species are different. For example, in the case of NV 
and PI spins, the zero-field splitting of the NV creates 
an energy mismatch. 

The swapping Hamiltonian can be reintroduced by 
inducing a Hartman-Hahn matching of the energies in 
the rotating frame under a continuous microwave irra- 
diation [31, 32]. By adjusting the Rabi frequency and 
the offset, the two spin species are brought into reso- 
nance and spin flip-flops (allowing polarization transfer) 
are now allowed, leading to a buildup of polarization. 
The environment spins can then be polarized efficiently 
by alternating periods during which the probe spin is po- 
larized and periods during which polarization exchange 
is driven by the microwave irradiation. 

The buildup of polarization can happen either via di- 
rect interaction between the probe spin and a spin in 
the environment, or indirectly via spin-diffusion [33, 34]. 
Since we are interested only in polarizing strongly cou- 
pled spins, the first process is dominant. Then, we can 
estimate the polarization time by the number of spins 
we want to polarize divided by their average coupling 
strength, T po i ~ n sc {T)/ (A) w n sc (T)T 2 * (where we 
used 1/T 2 * = yzT^fc to estimate the average coupling 
strength, an upper bound for T po i would be more gener- 
ally T pol < n sc (T)T/Ti). 

A different strategy to initialize the spin environment 
is measurement-based polarization [35] with cither feed- 
back or adaptive schemes. Precise measurement of the 
local magnetic field created by the spin environment at 
the sensor spin location effectively determines the envi- 
ronment spin state, with an increasing knowledge of the 
magnetic field shift corresponding to a reduced spin-state 
distribution and hence higher polarization. 

The polarization time will reduce the achievable sen- 
sitivity per root Hz, rj, since it increase the preparation 



time such that fewer measurement can be performed dur- 
ing a certain time interval. The exact sensitivity degrada- 
tion will depend on many factors, e.g. the depolarization 
(Ti) time of the environment spins, which determines 
how often the preparation step needs to be repeated. 

B. EAM with no polarization and phase error 

In the discussion so far we assumed that the external 
field to be measured was either static or oscillating in 
phase with the control sequence. For b(t) = b cos(27rf/r), 
we obtained the signal S = — sin$), where the phase 
is given by Eq. (3) for the EAM scheme or by $ = 2jbT/ir 
for the spin-echo scheme. However, if the field has a 
random phase (or cannot be synchronized perfectly with 
the pulse sequence) the signal averaged over many runs 
goes to zero, as (S) = |(1 - (sin$)) = if ($) = 0. 
Furthermore, even higher momenta of the signal, (S n ) 
are zero; thus it is not possible to infer information about 
the stochastic field by this method. 

A possible solution is to change the phase of the final 
pulse [36] (or equivalently, to introduce an additional, 
known phase accumulation during the free evolution). 
Then the signal becomes 

S=i[l + <cos($ + 0)>] (11) 

where $ is the phase difference between the initial and 
final pulse, $ is the phase due to the field to be mea- 
sured and we neglect any decay for simplicity. Since 
(cos($ + ■&)) = (cos$) cos the maximum signal is ob- 
tained for $ = 0, or by setting the phase of the initial 
and final pulse to be equal. 

The phase $ acquired in the modified EAM scheme 
(Fig. 2 with the last pulse along x) is different than that 
obtained in Eq. (3). In the limit of small fields, we obtain 
the signal 

S x = l- (2 + E fc [l + cos (^i) 2 ] sin (^) 2 ) 

«i-*(£) a [2 + N-], 

(12) 

where again we only summed over the "strongly coupled" 
environment spins. We note that the signal does not de- 
pend on the polarization of the environment spins (at 
least to first order in the polarization and to second or- 
der in the field b). Thus, even in the absence of any 
polarization it is possible to measure the external field 
(although not the sign of it). 

We compare the achievable sensitivity of the EAM and 
spin-echo method in the case where no polarization is 
present and the control sequence describe above is used. 
Optimizing the sensitivity with respect to the interroga- 
tion time r, we obtain the sensitivity for the spin echo 
sequence 

(13) 
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while for the EAM sequence we have 

Veam,x ~ , (14) 

CTV^y 1 + l n sc 

In the case of zero polarization (or of a signal with a ran- 
dom phase) it is no longer possible to obtain a quantum 
enhancement and have a scaling proportional to 1/ n sc by 
exploiting the spins in the bath. Indeed if there is no po- 
larization, no entanglement is created in the system, and 
no quantum enhancement of the sensitivity is expected 2 . 
Nevertheless, for favorable conditions of the spin envi- 
ronment (high quality Q and low ratio r) it might still 
be beneficial to use the EAM scheme instead of a sim- 
ple spin-echo magnetometry because it allows for an im- 
provement ~ ^/n sc by exploiting the unpolarized spins. 

VI. EXTENSION TO OTHER SPIN PROBES 

In the previous sections we presented a scheme that 
relied on the fact that for one of the eigenstates of the 
probe spin (|0)) the couplings to the environment spins 
was zero. It is possible to extend the EAM scheme to 
the case where the probe is a spin- 1/2, but only if the 
environment-probe spin couplings are all of the same 
sign. Such a situation could for example be realized by 
considering a single quantum dot in the nuclear spin en- 
vironment. The interaction between the central spin and 
the environment spins in this system is given by the con- 
tact interaction, whose strength depends mainly on the 
electronic spin wavefunction density and does not present 
the strong angular dependence of the dipolar interaction. 

To apply the EAM scheme with a spin-1/2 probe, we 
rotate the environment spins to be aligned along the y 
axis before applying the sequence shown in Fig. (2). For 
small fields, the additional phase acquired thanks to the 
environment spins is given by 




If all the couplings A& are positives, the environment 
spins contributions add constructively and it is always 
possible to find a time r s.t. there are n sc strongly 
coupled spins for which < AfcT < 47r, so that 

£ fc sin(Y) 3 «3^ S c 

More generally, it is also possible to use probes with 

higher spins, selecting two of their eigenstates as the lev- 
els of interest, by driving transitions on resonance with 
their energy difference. If the two eigenstates \a) and \b) 
are such their eigenvalues have different absolute values, 
\m a \ =/= \mb\ then we can apply the EAM sequence for any 



value of the couplings, as the phase enhancement will be 
^ ^ Sfe sm y ^" a ^~^ m ^ AfcT^j . Otherwise, one might use 

the modified scheme just presented in this section, if all 
the coupling constants are positive. 

As shown in this section, the scheme we introduced is 
quite flexible and can be applied to many different physi- 
cal systems, beyond the one we focused on in this paper. 
Besides spin systems, the same ideas could for example 
find an implementation based on trapped ions [2]. 



VII. DISCUSSION AND CONCLUSION 

In conclusion, we analyzed the EAM scheme intro- 
duced in Ref. [2] , which aims at enhancing the sensitivity 
of a single solid-state spin magnetic field sensor, by ex- 
ploiting the possibility to coherently control part of its 
spin environment. The environment spins act as sensitive 
probes of the external magnetic field, and their acquired 
phase is read out via the interaction with the sensor spin. 
Since the measurement scheme maintains roughly the 
same coherence times of spin-ccho-based magnetometry 
and the noise is still the shot-noise of a single qubit, we 
achieve a quasi-Heisenberg limited sensitivity enhance- 
ment. We analyzed in detail the sources of decoherence 
and confirmed with numerical simulations that the sen- 
sor coherence time under the EAM scheme is comparable 
to the Ti time under spin-echo, since the leading cause 
for decoherence has the same origin in the two cases. 
We further showed that dynamical decoupling schemes 
aimed at increasing the correlation time of the spin envi- 
ronment, by reducing the effects of intra-bath couplings, 
can be embedded in the measurement scheme and leads 
to longer coherence times and enhanced sensitivity. We 
extended the EAM scheme to the case where the envi- 
ronment spins are in a highly-mixed (zero-polarization) 
state, by appropriately modifying the detection sequence. 
This modified scheme achieve a classical scaling of the 
sensitivity, but can still be beneficial whenever the polar- 
ization methods we outlined cannot be applied or the AC 
field to be measured has a random phase. Our analysis 
finds that the sensitivity is determined by the "quality" 
of the environment, a parameter that takes into account 
a compromise between the number of strongly coupled 
environment spins with the reduced coherence time they 
entail. This result can be used to define the specifica- 
tions of engineered systems with controlled densities of 
spin impurities for optimal sensitivity. 
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